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Motivated by the ambiguous experimental data for the superconducting phase in silane (S 1 H 4 ), 
which may originate from the platinum hydride (PtH), we provide a theoretical study of the super¬ 
conducting state in the latter alloy. The quantitative estimates of the thermodynamics of PtH at 100 
GPa are given for a wide range of the Coulomb pseudopotential values (/z*) within the Eliashberg 
formalism. The obtained critical temperature value (Tc £ (12.94, 20.01} for fi* £ (0.05, 0.15}) agrees 
well with the experimental Tc for S 1 H 4 which may be ascribed to PtH. Moreover, the calculated 
characteristic thermodynamic ratios exceed the predictions of the Bardeen-Cooper-Schrieffer theory, 
implying occurrence of the strong-coupling and retardation effects in PtH. We note that our results 
can be of high relevance for the future studies on hydrides. 
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I. INTRODUCTION 

Recent experimental results on the superconducting 
phase in H 2 S hydride at high pressure [I], which report 
record high critical temperature value (Tc ~ 190 K), 
may become an another breakthrough in the research on 
the room-temperature superconductivity. In addition to 
the exceptional thermodynamic properties of the H 2 S al¬ 
loy 12], 0, 0, discussed material is considered to be 
a phonon-mediated superconductor. The latter fact is 
of great importance since electron-phonon superconduc¬ 
tors can be described within well-established Bardeen- 
Cooper-Schrieffer theory (BCS) 0, (6( and its derivatives 
such as the strong-coupling Eliashberg formalism .7], 

In this spirit, hydride materials, whose concept stems 
from the milestone works by Ashcroft 0 , m , are cur¬ 
rently one of the most promising candidates for the fu¬ 
ture high-temperature superconductors nni in compar¬ 
ison with e.g. cuprates mi, m- Among the family 
of hydrogen-based materials, the high-pressure supercon¬ 
ducting phase has been experimentally observed only in 
the aforementioned H 2 S alloy and likely in the platinum 
hydride material (PtH). Please note however that these 
results are not fully confirmed yet. 

From the scientific point of view, the irrefutable testi¬ 
monies for the existence of the superconducting phase in 
the PtH and H 2 S alloys are equally important for push¬ 
ing forward research in the discussed domain. However, 
in the case of PtH the already available results are not 
so clear and direct like in the case of the H 2 S hydride. 

In particular, in 2008 Eremets et al. conducted the 
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high-pressure experiment for silane (SiH 4 ) and stated 
that the induction of the superconducting state in this 
material is possible [13] . However, further investigations 
of the results obtained by Eremets et al. questioned to 
some extent the outcome of the discussed work. Degt¬ 
yareva et al. claimed in their paper m that, in the 
discussed experiment, silane can decompose and released 
hydrogen may react with the surrounding metals, in par¬ 
ticular with platinum. This statement was motivated by 
the XRD data which matches PtH instead of the SiH 4 al¬ 
loy, suggesting that the observed superconducting state 
may originate from the first mentioned material. Further 
theoretical studies confirmed that such scenario is highly 
possible m, m, m, cu However, all mentioned the¬ 
oretical papers presented approximate estimations of the 
critical temperature by using the McMillan formula [19], 
and did not provide the discussion of other crucial ther¬ 
modynamic properties. In this context, the discussed 
problem of the existence of the superconducting state in 
PtH is still an open and emerging question. 

In order to address raised in the previous paragraph is¬ 
sues, we present complementary and comprehensive the¬ 
oretical analysis of all of the most important thermody¬ 
namic properties of the PtH superconductor at the pres¬ 
sure of p = 100 GPa. This pressure value is chosen in 
order to match the conditions at which the highest value 
of the critical temperature has been obtained in the ex¬ 
periment of Eremets et al. [T5], In particular, our calcu¬ 
lations are conducted within the Eliashberg theory due 
to the relatively high electron-phonon constant for the 
PtH hydride (A = 0.84 as predicted by Kim et al. |T8j). 

The present paper is organized as follows, in the first 
section we briefly elaborate on the technical aspects of the 
Eliashberg formalism and applied numerical procedures. 
Next we present our numerical results obtained on the 
imaginary axis as well as in the mixed representation 
of the Eliashberg equations. Our paper ends with some 
pertinent conclusions. 
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II. THEORETICAL MODEL 

The theoretical analysis presented in this paper based 
on the isotropic Eliashberg equations which are solved, 
at the imaginary axis and in the mixed representation, 
by using the iterative method given in [50] and [21] . The 
isotropic character of the Fermi surface is imposed by 
the Eliashberg spectral function (a 2 E(0)), which mod¬ 
els the electron-phonon interactions and was taken from 
m- We note, that in order to assure the consistency of 
our calculations with the experimental conditions in m , 
the chosen a 2 F{ fi) function corresponds to the pressure 
value of 100 GPa and the hep structure of the space group 
P6 3 . 

Due to the uncertain experimental data we consider 
three distinct values of the Coulomb pseudopotential (p*) 
which models the depairing electron correlations and al¬ 
lows us to sample the superconducting phase of PtH. In 
particular, p* = (0.05,0.1,0.15) and match the general 
estimations suggested by Ashcroft for the hydride super¬ 
conductors in |8j. Furthermore, the wide range of the p* 
values allows us to appropriately examine convergence of 
our results and the experimental data of Eremets et al. 
na. as well as investigate the influence of the p* value 
on the superconducting properties of PtH. 

The precision of our numerical analysis is controlled 
by assuming the sufficiently high limit for the Matsub- 
ara frequencies: w m = ^(2m — 1), where j3 = 1/fcgT 
is the so-called inverse temperature, and kg states for 
the Boltzmann constant. Our calculations are conducted 
then for the 2201 Matsubara frequencies, allowing us to 
estimate quantitatively the superconducting properties 
of the PtH for temperatures T > T 0 = 3 K. Finally, the 
cut-off frequency is set to be to c = 10H max , where f2 max is 
the maximum phonon frequency equals to 193.25 meV. 

III. NUMERICAL RESULTS AND DISCUSSION 

In the first step we analyze the temperature depen- 
dance of the order parameter (A m ) by conducting our 
numerical calculations on the imaginary axis. The ob¬ 
tained results are presented, for all considered p* values 
and selected temperatures, in Figs, [l] (A-C) as a func¬ 
tion of m. The physically reasonable gradual decrease 
of the A m value with increasing p* and T can be ob¬ 
served. Moreover, together with increasing m we notice 
saturation of the obtained solutions in all cases. It is a 
testimonial for the precision of our calculations and the 
quantitative character of the obtained results. 

The value of the critical temperature is determined 
on the basis of the condition: A m= i(Tc<) = 0, where 
A m —i denotes the maximum value of the order parame¬ 
ter. Supplementing functional dependance of the A m= i 
function on the temperature for the selected values of p* 
is depicted in Fig. [I] (D). In particular, calculated Tc 
amounts 20.01 K, 15.78 K, and 12.94 K for p* equals to 
0.05, 0.1, 0.15, respectively. 




FIG. 1: (A)-(C) The order parameter (A m ) on the imaginary 
axis as a function of m for selected values of the temperature 
( T) and Coulomb pseudopotential ( p *). (D) The maximum 
value of the order parameter (A m= i) on the imaginary axis 
as a function of the temperature for the selected values of the 
Coulomb pseudopotential. 


We note that the experimental value of Tc obtained 
by Eremets et al. [T3] is ~ 17 K for p ss 100 GPa. In 
what follows, it matches our predictions falling in to the 
physically feasible region of p* £ (0.05,0.1). 

The results obtained on the imaginary axis and pre¬ 
sented in Figs. |T] (A-D) provide also information on 
other measurable observable, namely the effective elec¬ 
tron mass (?n*). Specifically, the value of the to* can be 
determined by using the following relation: 

to* ~ Z m=1 m e , (1) 

where Z m= \ is the maximum value of the wave function 
renormalization factor and m e denotes the band electron 
mass. 

The wave function renormalization factor as well as 
its maximum value are presented in Figs. [2] (A-D) as a 
functions of the temperature for the selected values of 
p*. By using Eq. ([l]) we estimate the to* value for the 
PtH alloy to be 1.84 m e , which is relatively high as for 
the phonon-mediated superconductors. 

The determined dependance of the Z m= \ function on 
the temperature allows us next to calculate the normal¬ 
ized free energy difference between the superconducting 
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FIG. 2: (A)-(C) The wave function renormalization factor on 
the imaginary axis (Z m ) as a function of m for the selected 
values of temperature (T) and Coulomb pseudopotential (p*). 
(D) The maximum value of the wave function renormaliza¬ 
tion factor on the imaginary axis {Z m = 1 ) as a function of the 
temperature for the selected values of the Coulomb pseudopo¬ 
tential. 


and normal state as: 
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In Eq. ([2]), p( 0) states for the electron density of states 
at the Fermi level, whereas and Z^ are the wave 
function renormalization factors for the superconducting 
(5) and normal (N) state, respectively. 

The AF/p(0) as a function of temperature for three 
considered p* values is presented in the lower panels of 
Figs. [3] (A-C). We observe that in all cases AF/p(0) 
takes negative values indicating thermodynamic stabil¬ 
ity of the superconducting phase, where togheter with 
decreasing value of p* the stability increases. 

Figs. [3] (A-C) present also corresponding functional 
behavior of the normalized thermodynamic critical field 
Hc/\Jp{ 0) on temperature, obtained on the basis of the 
following expression: 


H c 

VpW) 


v/- 8 7T [AF/p(0)]. 


( 3 ) 


Like in the case of the AF/p(0) function, the normal¬ 
ized thermodynamic critical field strongly depends on the 
assumed p* value, and decrease together with increasing 


p*. The F[c /\Jp (0) function and the specific heat for 
the normal state (C N ) given by: 

° N m = 1 (4) 

k BP ( 0) 0' 1 J 


where 7 = | 7 r 2 (1 + A) is the Sommerfeld constant, can 
be next used to to estimate the first characteristic ther¬ 
modynamic ratio, familiar in the BCS theory. In partic¬ 
ular, the ratio of interest for the critical magnetic field is 
defined as: 


= T c C n ( T c ) 

Hi ( 0 ) 


( 5 ) 


and takes the constant value of R|j cs =0.168 within the 
BCS theory [5], 6 ]. However, for the PtH alloy at 100 
GPa our calculations predict that Rh £ (0.155,0.167) 
for p* £ (0.05,0.15), showing deviation from the BCS 
predictions. 




t (K) 


T (K) 


FIG. 3: (A-C) The normalized free energy difference 

(A F/p (0)) (lower panel) and the thermodynamic critical field 
( Hcl\JP (0)) (upper panel) as a function of the temperature 
(T) for the selected values of the Coulomb pseudopotential 
(M*)- 


Relations given in Eqs. © and ([9]), can be next used 
to address the normalized difference between the specific 
heat of the superconducting (C 5 ) and normal (C^) state, 
as given by: 

A C(T) _ 1 d 2 [AF/p (0)] 

k B cp{ 0) 0 d (k B T) 2 ’ 1 

wheie AC = C s — C N . In what follows, the specific 
heat of the superconducting state can be determined. In 
Figs. [3] (A-C) we present obtained functional behavior 
of the C s on the temperature for the three considered 
p* values. For comparison and systematical purposes, 
the C N function is also presented. It can be observed, 
that C s function rises faster then C N as temperature 
grows, presenting characteristic jump at Tc■ Moreover, 
the normalized specific heat of the superconducting state 
strongly depends on the p* value, decreasing when the 
p* grows. 
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For the results given in Figs. [4j the second dimen¬ 
sionless characteristic thermodynamic ratio can be cal¬ 
culated, by using: 


= A C(T c ) 
C N ( T c )' 


(7) 


Similarly like in the case of the Rh ratio, the obtained 
Rc values exceed the BCS limit. In particular, Rc G 
(1.88,1.82) for p* G (0.05,0.15), whereas R% cs = 1.43. 




T (K) T (K) T (K) 


FIG. 4: The normalized specific heat of the superconducting 
{C s /ksp{ 0)) and normal [C N /ksp{ 0)) state as a function 
of the temperature (T) for the selected values of Coulomb 
pseudpotential (p*). 


The last remaining characteristic ratio, which corre¬ 
sponds to the energy gap value at the Fermi level (2A(T), 
where A (T) is the the order parameter on the real axis), 
can be calculated by solving the Eliashberg equations 
in the mixed representation (please see [20] and [21] for 
more details). Such analysis is required in order to pro¬ 
vide the quantitative values of the A (T) function, given 
by the following expression: 


A (T) = Re [A (w = A (T), T)\. (8) 

In Fig. [5| we summarize our calculations on the real 
axis by presenting the results of the total normalized den¬ 
sity of states (NDOS(w)) for T = Tq = 3 K and the se¬ 
lected values of Coulomb pseudopotenital. We note, that 
the NDOS(w) function is related to the A (T) solutions 
as: 



<o (meV) 


FIG. 5: The normalized total electronic density of states 
(NDOS) at To = T = 3 K for the selected values of Coulomb 
pseudopotential (p*). 


of p *, we can clearly observe the decrease of the supercon¬ 
ducting gap towards metallic character of the discussed 
alloy. This outcome is in qualitative agreement with the 
results presented previously in Figs. S(A-D). 

On the basis of the obtained results we first estimate 
the value of the the zero temperature energy gap at the 
Fermi level (2A(0) = 2A (To)). In the case of the PtH 
alloy at the pressure of 100 GPa the calculated values are: 
2A(0) G (6.87,5.31,4.28) meV for p* G (0.05,0.1,0.15). 
In the context of the BCS theory, the 2A(0) parameter 
can be next used to determine the thermodynamic ratio 
for the zero temperature energy gap at the Fermi level 
(Ra)- In particular: 


Ra 


2A(0) 

k B T c ' 


( 10 ) 


Again, the BCS theory predicts constant value of the 
Ra ratio equals to 3.53, whereas our results for the 
PtH material under discussed conditions give Ra G 
(3.98,3.84) for p* G (0.05,0.15). In what follows, sim¬ 
ilarly like in the case of the R B and Rc ratios, the Ra 
differ from the BCS suggestions. 


NDOS (w) 


DOS s (w) _ 
DOS n (w) “ 


| u) — ?T| 

\A W - * r ) 2 - A 2 M 


(9) 

where DOSs (w) and DOSn (w) are respectively the den¬ 
sity of states of the superconducting and normal state, 
whereas T states for the pair breaking parameter and 
equals to 0.15 meV. 

As already mentioned, results presented in Fig. [5] are 
depicted for T = T 0 due to the fact that only the maxi¬ 
mum value of order parameter is relevant for our discus¬ 
sion. By comparing plots for different considered values 


IV. SUMMARY 

In the present paper we address the vague experimen¬ 
tal data for the high-pressure silane superconductor given 
by Erements et al. in |13] , and the consequences of these 
results in the terms of the platinum hydride alloy. In par¬ 
ticular, we perform calculations of the thermodynamic 
properties of the PtH hydride and reveal several crucial 
facts concerning possibility of the existence and proper¬ 
ties of the superconducting phase in this material. In or¬ 
der to assure the creditability of our resutls, calculations 
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are preformed at the pressure of 100 GPa and within the 
hep structure (with the P 63 space group), which match 
exactly the conditions of the experiment conducted in 

m- ' 

In the first place, the calculated values of the critical 
temperature for three different p* present good agree¬ 
ment with the experimental data. In particular, it is 
shown that the experimental value of Tc can be obtained 
for the p* equal to ~ 0.1. This is physically relevant 
value, which additionally agrees with the predictions of 
Ashcroft given for the hydrogen-dense materials in | 8 j. In 
what follows, our results indicate that experimental data 
obtained in m, with a high probability, can be ascribed 
to the PtH alloy rather then S 1 H 4 material. This outcome 
is in agreement with the previous qualitative results for 
Tc presented in [T5], 1(1 . [17], [TSJ. 

Our further calculations, allows us to determine other 
thermodynamic properties of the PtH at 100 GPa such 
as the thermodynamic critical field, the specific heat for 
the superconducting and normal state, the value of the 
energy gap at the Fermi level, and others. 

It is convenient to summarize these results by esti¬ 
mating the values of the characteristic thermodynamic 


ratios. We show that all calculated ratios exceed pre¬ 
dictions given by the BCS theory. In particular Rh £ 
(0.155,0.167), R c £ (1.88,1.82), and R A £ (3.98,3.84), 
which correspond to the range of p* £ (0.05,0.15). The 
differences between our estimations and the Bardeen- 
Cooper-Schrieffer theory indicate a significant role of the 
retardation and strong-coupling effects in PtH at the 
pressure of 100 GPa. In what follows, the analysis of 
the discussed material cannot be conducted within the 
BCS theory, reinforcing the relevancy of our results and 
chosen theoretical approach. We note also, that the in¬ 
crease of the value of the Coulomb pseudopotential re¬ 
duces mentioned discrepancies, which is the physically 
plausible effect. 
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